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Abstract 

We use the Cole-Hopf transformation and the Laplace method for the heat equa- 
tion to justify the numerical results on enstrophy growth in the viscous Burgers 
equation on the unit circle. We show that the maximum enstrophy achieved in the 
time evolution is scaled as £ 3 / 2 , where £ is the large initial enstrophy, whereas the 
time needed for reaching the maximal enstrophy is scaled as £~ 1 / 2 . These bounds 
are sharp for sufficiently smooth initial conditions. 



Existence and regularity of solutions of the three-dimensional Navier-Stokes equations 
is a challenging problem that attracted recently many researchers [7J HO]. One possibility 
to continue local solutions globally in time is to control enstrophy of the Navier-Stokes 
equations during the time evolution. Enstrophy is an integral quantity in space, which may 
diverge in a finite evolution time, indicating singularities in the Navier-Stokes equations. 
To control the growth of enstrophy, Lu and Doering [TT] studied bounds on the instanta- 
neous growth rate of enstrophy and showed numerically that these bounds are sharp in the 
limit of large enstrophy. However, these bounds on the instantaneous growth rate do not 
imply that enstrophy blows up in a finite time, because solutions of the Cauchy problem 
associated with the Navier-Stokes equation may deviate away from the maximizers of the 
bounds even if the initial data are close to the maximizers. 

To deal with this problem, Ayala and Protas [2] looked at a toy model, the one- 
dimensional viscous Burgers equation. Their numerical results indicated that the bounds 
on the instantaneous growth rate are not sharp when they are integrated on a finite time 
interval in conjugation with the energy dissipation and the Poincare inequality. Limited 
accuracy of numerical results did not allow them to conclude if better estimates on the 
enstrophy growth over a finite time interval can be justified within this context. 
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To address this question, we consider the Cauchy problem for the one-dimensional 
viscous Burgers equation [3J, 



u t + 2uu x = u xx , x g T, t e R+, , . 

u\ t=0 = u , x G T, 

where T = [— |, |] is the unit circle equipped with the periodic boundary conditions for 
the real-valued function u. Local well-posedness of the initial-value problem (TjQ) holds for 
«o £ H^ er {T) with s > — | [6]. Global existence holds in Hp ev (T) for any integer s > 0, 
and we consider here global solutions of the viscous Burgers equation flTJ in H^ eT (T). For 
simplicity, we will assume that uq is odd in x, which implies that u(—x, t) = —u(x, t) holds 
for all t > and i6l 

All Ljp er norms with even p are monotonically decaying in the time evolution of the 
viscous Burgers equation (JTJ. In particular, the energy dissipation equation follows from 
([I]) after integration by parts: 

If, dKiu) f . 

K[u) — — j u ax =>■ — - — = / u(u xx — 2uu x )dx = —2E(u), (2) 



2 Jf dt 

where E(u) = \ J T u 2 x dx is a positive definite enstrophy. For a smooth solution u G 
C(R + , Hp er (T)), the enstrophy changes according to the equation: 



dE{u) 
dt 



j u x (u xxx - 2uu xx - 2u 2 x )dx = - (u 2 xx + ul)dx = R(u), 



(3) 



where R{u) is the rate of change of E(u). We can see from that R(u) is a sum of negative 
definite quadratic part and a sign- indefinite cubic part. The quadratic part corresponds 
to the diffusion term of the viscous Burgers equation and the cubic part corresponds to 
the nonlinear advection term. It is the latter term that may lead to the enstrophy growth 
during the initial time evolution. 

Lu and Doering jTTJ showed that the rate of change R(u) in (J3j can be estimated by 

R(u) < \E"\u\ (4) 

and illustrated numerically that the growth R{u) = 0(S 5 ^ 3 ) is achieved in the limit of 
large S := E(u). If bound (jl]) is sharp on the time interval [0,T] for some T > 0, then 
integration of the enstrophy equation ([3]) with the use of the energy dissipation equation 
implies 

EV\u(T)) - E^(u ) <\j^ E(u(t))dt = ~ [K(u ) - K(u(T))\ . (5) 
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Using the Poincare inequality for periodic functions with zero mean, 

K(u ) < (6) 
and neglecting K{u{T)) in (jSJ), we can obtain 

E(u(T))< ^/3 + _i_^ 3 5 ^ :=jE ; M . (7) 

Ayala and Protas [2] showed numerically that the integral bound (jTJ) is not sharp even in 
the limit of large S. Instead, they obtained the following numerical result: 

T* = 0(S^ 2 ), E(u(T*)) = 0(£V 2 ), K(u(T*)) = 0(£), (8) 

where is the value of time t, at which E(u(t)) is maximal. They also wanted to show 
that K(u ) — K{u{T*)) = 0{S 1 ^ 2 ), so that the full integral bound ^ could be useful but 
numerical approximations of this quantity suffered from large errrors: 

K(u Q ) - K{u{T*)) = O(£°- 68±0 - 25 ). (9) 

In the previous work [13], we used dynamical system methods to study analytically the 
enstrophy growth in the viscous Burgers equation. Our technique based on the self-similar 
transformation and analysis of asymptotic stability of viscous shocks in an unbounded 
domain did not rely on the remarkable properties of the viscous Burgers equation such as 
the Cole-Hopf transformation [5] [H] of equation ([T]) to the linear heat equation. On the 
other hand, a weaker version of bounds OH]) modified by logarithmic factors was justified 
as a result of this approach. No estimate on K(uq) — K(u(T*)) has been obtained. 

In this work, we shall rely on the Cole-Hopf transformation and use the Laplace method 
for the heat equation. The Laplace method is typically used to recover solutions of the 
inviscid Burgers equations from solutions of the viscous Burgers equation in the limit of 
vanishing viscosity (see, e.g., [15] Chapter 2], [TJ Example 6.5.2], or [12] Section 3.6]). 
Applications of this method to statistical properties of the Burgers turbulence can be 
found in [8]. 

Note that the limit of vanishing viscosity corresponds to the limit of large enstrophy 
in the context of our work. We shall implement the Laplace method to justify numerical 
results OH]) and to estimate K{uq) — K{u(T*)) as £ — > oo. Our main result is the following 
theorem. 

Theorem 1 Consider the initial-value problem (QP with initial data uq{x) = kf(x), where 
f G Cp Cr (T) is an arbitrary odd function such that f"(x) > for all x G [0, |] . Consider 
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the limit k — > oo and denote the initial enstrophy by £ = E(u ) = 0(k 2 ). There exists 
> such that the enstrophy E(u) achieves its maximum at = u(-,T*) with 

T, = 0(£^ 2 ), E(u*) = 0(£ 3/2 ), K(u*) = 0(£), (10) 

and 

K{u )-K{u*) = O(£), (11) 
where all bounds are sharp as £ — > oo. 

Because / G Cp er (T) is an odd function with f"(x) > for all x G [0, |], it implies 
necessarily that /(0) = / (~) = 0, f(x) < for all x G (0, |), and f'(y) is a monotonically 
increasing function from f'(0) < to /' (|) > with a unique zero at such that 
/'(#*) = 0. We will actually show that 

r * = 2fc|/(s,)l ' ^) = ^ 3 !/(^)| 3 + ^ 2 ), as fc^oo, (12) 

and 

K(u )-K(u*) = k 2 ^ * f( y )dy~^xj 2 (x^ + 0{k), as k ^ oo, (13) 

where the leading-order terms are all nonzero. Sharp bounds (fT2l -(fT3l) rule out the hope 
of using the integral bound (jSJ) that follows from the instantaneous estimate (BJ and the 
balance equations (j2J) and (j3J). 

As an example, we can consider f(x) = — 2n sin(27rx), which saturates the Poincare 
inequality and satisfies the conditions of Theorem [IJ Then, x* — ^, where /'(^*) = 0, 
and the sharp bounds (IT2"1) - ([T3"1) yield the explicit expressions 

1 7T 2 

7; = ——, £«) = 4^ 3 A; 3 + e>(fc 2 ), if (%) - = —A; 2 + as k -> oo. 

lo7rA; o 

A generalization of Theorem [T] can be developed for any / G Cp Cr (T) with finitely many 
changes in the sign of /" on [0, |] by the price of lengthier technical computations. On the 
other hand, it is not clear if the bounds (fT0|) and ffTTj) remain sharp for initial conditions 
with limited regularity. 

In the remainder of this paper, we shall prove Theorem [TJ Using the Cole-Hopf trans- 
formation [5J |H] , 

d 

u(x, t) = — — \ogip(x, t), mx, t) > for all (x,t), (14) 
ox 

we rewrite the Cauchy problem ([1]) in the equivalent form 

4>t = ipxx, x G T, t G R + , . lg v 

*0|t=o = </>o, ac e T, 
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where 

^ (x) = e -J?«oW* = e - fcF W, := / /(a)ds. 



If / G Cp er (T) is odd in x, then ^ G CperC^) * s even in x. 

Although the Cauchy problem (|15|) is posed on the periodic domain, we can still con- 
struct the solution as a convolution of the initial data ipo with the heat kernel Gf : M — > M + 
defined on the entire axis, 

1 x 2 

G t (x) = -=e"tf, xeR, t G R + . (16) 
This is justified by the generalized Young inequality, 

\\G t *i>o lU^M) < II L p(m) || V'oIIjwcr), l<P,g,^<oo: - + - = 1 + -, 

where * denotes the convolution integral. If G -^ r (^)> then G 4 * -0 O G L^ r (M) for all 
t G R + . Therefore, we write the solution to the Cauchy problem (Tl5|) in the explicit form, 

i;(x,t) = ^= e- kF ^-^dy. (17) 

Let us define the parametrization of the time variable by t — ^r-, where a G M+ is a 
new parameter. The solution of the viscous Burgers equation can now be written in the 
explicit form, 



u(x,t) — — ^- log I X:a (k), I x , a (k) := J 



oo 

p -ktj) x , a {y) 



dy, (18) 



where (f) x ,a(y) '■= F(y) + ~a(x — y) 2 . Integral I x , a (k) can be studied in the limit k — > oo 
by means of the Laplace method (Section 3.4 in |12j). We recall the main result of the 
Laplace method. 

Proposition 1 For any — oo < a < b < oo, assume that (p G C 4 (a, b) has a global 
minimum at c G (a, b) such that <fi'(c) = and <p"{c) > 0. Then, for any 9 G C 2 (a,b) n 
L 1 (a, b), we have the following asymptotic expansion: 



b / ~ \ 1/2 



J(A;) := / 9{y)e- k *^dy = ( ) %)e-^ c ) 



as /c — )■ oo. (19) 
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Here and in what follows, we use the following notations. Let X be a Banach space. 
We write A = Ox{k p ) as k — > oo if there exist constants C± such that < C_ < C + < oo 
and C„£ p < \\A\\ X < C + £ p . If X = R, we write A = 0{£ p ). 

Applying Proposition (TJ we obtain the following. 

Lemma 1 Under assumptions of Theorem^ the following expansion 

u(x,t) = kf(s Xj a) + 0L~ r (T)(l) as k -> oo, (20) 
holds for all t G [0,£o) ; where t := 2 k\f'(o)\ an d s x , a is the unique root of a(s — x) + f (s) = 0. 



Proof. We compute 

^,a(y) = f(y) + a (y - <f>xM = f'(v) + a - 

For every a > max ye f (—f'(y)) = |/'(0)| and every x G T, there exists exactly one root 
of f(y) + a(y — x) = 0. Let us denote this root by s x>a . Conditions of Proposition [1] are 
satisfied, so that 



Ix,a(k) 



2n 



k(f'(s x>a ) + a) 



1/2 



&0a;,a(Sx,a) 



as k — > oo 



and 



d x Ix,a(k) = ak / (y-x)e- k *"^dy 



ak^Sx^a X) 



2tt 



k(f'(sx,a) + a) 



1/2 



1 + 0| i 



as — > oo. 



Using these expansions, the representation (fT8l) . and the equation a(x — s x , a ) = f(s X)a ), we 
obtain (|2Q). □ 

Remark 1 TTie result of Lemma [J\ is known in the limit of vanishing viscosity, when 
a smooth solution of the viscous Burgers equation is shown to converge to the classical 
solution of the inviscid Burger equation before a shock is formed JTS[ Chapter 2j. 



Using the representation 



u x (x, t) = u 2 (x, t) — 



Ix,a(k) 



(21) 



and the Laplace method for d^I Xja (k), the result of Lemma [T] can be extended to show for 
all t G [0,t ) that 

u x (x,t) = k a l^A + C> X oc r(T) (l) as A;^oc. (22) 
As a result, the energy and enstrophy are expanded as k — > oo as follows: 

= Ifc 2 / f(s x , a )dx + 0(k) 

and 

Because /(0) = / (|) = 0, we have s Xja = x for x = and x — \. Moreover, the map 
T 3 x t— > s X}(i G T is one-to-one and onto. This argument gives the energy conservation at 
the leading order for all t G [0,to) : 

K(u(t)) = Jj\y)[l + -J'{y)^dy + 0{k) 

= \k 2 J f(y)dy+0(k) 
= K(u ) + O(k). 

On the other hand, the enstrophy grows initially but remains within the 0(k 2 ) order for 
all t G [0,t ): 

E{u{t)) = \k* f^§fdy + 0(k). 
2 J T a + f\y) 

At t — to, that is, for a = max y£ j(—f'(y)) = |/'(0)|, a local (pitchfork) bifurcation 
occurs among the roots of f(s) + as = near s — 0. For t > t , /'(0) + a < 0, so that for 
all x near 0, a local maximum s x>a of exists near 0, whereas two local minima s^ a and 
s*, a of 0*,o exist for s+ a G (0, ~) and s~ a G (-|,0). 

In what follows, we consider the values of x G [0, |] and use the odd symmetry of 
u(x, t) for x G [— |,0l. Note that 0o,a(s^ a ) = 0o,a(% a ) due to the pitchfork symmetry 
When x is increased in (0, |), then 0cc,a( s J a ) < 0a;,a( s x tt ) an d the difference is growing with 
the values of x. Moreover, there is x G (0, |) such that the roots s x ^ a and s~ a coalesce at 
x = x and disappear as a result of the saddle-node bifurcation. To prove these claims, we 
denote 

¥x,a = </v(s~J - 4>xA S ta) = 2 fl ( S x,a ~ 4,a)(4,a + ~ 2x) - / /(s)dfi 

" J s x,a 



7 



and 

7'«J + ^ 1/2 

Xx,a 



f'{s- a ) + a 



We have the following. 



Lemma 2 For every t > t , where t := 2 k\f'(o) \ > ^ere is x G (0, |) such that for all 
x G [0, Xo), three roots of a(s — x) + f(s) = exists in the following order: 

1 ,1 

~2 < S *,a < S x,a < < S J >0 < -. 

Moreover, (p x>a and Xx,a are C 1 monotonically increasing functions of x with t^ ,a = 0, 
Xo,a = 1, and Xx,a — > +oo as x — > xq. The point xq marks the saddle-node bifurcation 
among the roots of a(s — x) + f(s) = such that s~ oa = s xo ^ a and a + f'(s X(ua ) = 0. 

Proof. The presence of three roots Sq a = —SQ a and So,a — follows for x = because 
a + /'(0) < for t > t and / is an odd function. By continuity, three roots persist 
and, by the implicit function theorem, the three roots are C 1 functions of x as long as 
a + f'{s xa ) > 0. We now compute 

d x (p x , a = a(s+ a - s"J > 0, 

and 



9xXx,a Xx,a 



>o, 



(« + /'(4J) 2 (a + ffe)) 2 ] 

where the last inequality is due to f"(x) > for x G [0, |] and the fact that /" is an odd 
function. In addition, we note that 

1 _ f4,a 

¥0,a = 2 a ( S 0,a - 4,a)(4,a + S 0,a) ~ j _ f( S ) d S = 



and 



, 7K) + a V /2 _ 1 



because s = — and / is odd. The statement of the lemma is proved. □ 

Remark 2 When a is reduced further, additional saddle-node bifurcations occur among 
the roots of f(s) + as = due to periodicity of f . For all x near 0, these bifurcations 
give rise to new maxima and minima of <f) X)a outside of T and the values of <f> x>a at new 
local minima are larger than the values of 4> x ,a(4,a)- Therefore, we can simply neglect the 
presence of these additional bifurcations in the applications of the Laplace method for all 
ae(0,|/'(0)|). 
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We shall apply Proposition [T] to obtain the following. 

Lemma 3 Let s^ a and x be described in Lemma [B Under assumptions of Theorem [7J 
for every t > t , where t := 2 h\f(o)\ > there is x\ G (0, x ) such that for all x G [0, X\], 

i ,\ , f( S x,a) + Xx,af{s x a )e lfix - a 

u{x,t) = k — — V L ^ {Q . Xl) {l) as k -> oo (23) 

whereas for all x G [xi, |] , 

M(x,t) = A:/«a) + LO c( !Blj ^(l) as k^oo. (24) 



Proof. We split the integral in the explicit solution ( |T8l) into two parts 



k**,a(v)d y+ / e~ k ^ a{y) dy 



and apply Laplace method of Proposition [T] separately for each integral resulting in the 
following 



2tt 



1/2 



+ 



Hf'(s- a ) + a) 
2tt 

fc(/'«a)+0) 



e 

1/2 



J + C "i 



1 + 0| i 



as k ^ oo 



and 

d x I x , a (k) = ak(s~ a - x) 
+ak(s^ a — x 



2n 



1/2 



*(/'(**,■) + ^ 

2tt 



1 + ol i 



1/2 



1 + ol i 



as — > oo. 



For x = 0, both leading-order terms in these expansions have equal magnitude resulting 
in (|2"3"|) . whereas for any fixed x > 0, the first integral is exponentially small compared to 
the second integral resulting in (1241) . □ 



Remark 3 The two different expansions |Hj) and (Efy are asymptotically equivalent for 
x = x\ (and near x = x\) because (p X}0 > for any x > and, therefore, the term involving 
e -kip x , a is exponentially small as k — >■ oo for any fixed x > 0. 
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Using the representation fl2Tl) and the Laplace method for d^I x>a (k), the result of Lemma 
[3] can be extended to prove for all t > t and all x G [0, xi]: 



u x (x,t) = -£r 



.Xx,a{f( S ta) - f( S x,a)) 



2p kipx^a 



+ Oi,ao( 0tXl \(k) as k — > oo, (25) 



(1 + x^e"^) 2 
whereas for all x G [xi, |] , 

u x (x,t) = (D LOB fi\(k) as fc ->■ oo. 
We shall now compute the leading order of the enstrophy for t > £ : 

E(u(t)) = / u 2 x (x,t)dx + 0(k 2 ), as fc -» oo. 



(26) 



(27) 



Because </? XiQ is monotonically growing for all x G [0, X\], we can use another version of the 
Laplace method (Section 3.3 in [T2]). 

Proposition 2 For any c > 0, assume that cj) G C 2 (0,c) is monotonically growing such 
that 4>'(0) > 0. Then, for any 9 G C 1 (0,c) fl L 1 (0,c), we have the following asymptotic 
expansion: 



Uk) := [ C 6(y)e- k ^dy = -M_ e -^(°) 

1 ; 7o ^(o 



as — >■ oo. 



(2? 



We have computed previously: ip ,a = 0, Xo,a = 1, and 

d x <p ,a = a(sja - So, a ) = 2twja = - 2 f(4,a) > °- 

Therefore, a straightforward application of Proposition [2] yields 

S(tt(0) = ^* S |/(«J a )| 3 + O(* a ), as fc^oo. (29) 

For a = |/'(0)|, 4 a = 0, and the £>(£; 3 ) term of (P2J vanishes because /(0) = 0. When a is 
decreased from |/'(0)| to 0, the value of s£ a grows from to | and it passes the value X*, 
where f'(x*) = and |/(sj" a )| is maximal. The corresponding value of a* is found from the 

equation = x # , or explicitly, a* = . This argument completes the proof of the 

first two bounds (HDD of Theorem Q] with T* = ^ = afclfa] ( reca11 that k = °(£ 1/2 ))- 

Remark 4 An application of the Laplace method to the values of t near t that depends 
on k is much more delicate and has been explored in the pioneer paper ^ (see [Lflj for 
recent development). We are very fortunate here that the main result of Theorem^ can be 
proven without knowing the behavior of the solution near t = t . 
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To complete the proof of Theorem [Q we need to compute the energy K(u(t)) for all 
t > t . By Lemma El we represent for all i£T: 

u(x,t) = kf(s+ a ) + u(x,t), 

where u is found from the asymptotic expansions ( 1231) and ( |24f) . By Proposition [2} we have 

yl/2 

K{u{t)) = k 2 f 2 (st a )dx + 0(k) as fc ->• oc. 
Jo 

For all a < |/'(0)|, we have s,j~ a > 0. Nevertheless, we still have s+ a = x for a; = |. The 
map [0, | 9i4 s+ € [ s o,ai §] i s one-to-one and onto. As a result, for all £ > i ; 

K(u(t)) = k 2 f 1/2 f(y)(l + \f(y)jdy + 0(k) 

^ ny)dy + ^\f(s+ a )A+0(k). 

0,a J 

This argument completes the proof of the third bound fflUj) of Theorem [TJ To prove the 
bound ( fTTj) . we need to show that the 0(k 2 ) term in the expansion for K(u(t)) is different 
from the one for K(uq) (it can only be smaller), or equivalently, that 

_i_ 

f{y)dy > ^\f(4, a )\ 3 = lsU 2 «a)- (30) 

To prove ( )30|) . we define two functions F, H : [0, x*] — >• R by 

G(x) := f f{y)dy - \xf\x), H(x) = f(x) - xf(x). 
Jo 

By the conditions on /, it is clear that G G C 3 ([0, x*]) and H G C 2 ([0, £*]). Furthermore, 
G(0) = 0, H(0) = 0, and 

G'(x) = \f{x)H{x), H\x) = -xf"(x). 

Because f"(x) > for all x G [0, 5], H is a monotonically decreasing function from 
-ff(O) = to if (a;) < for all x G (0,x*]. Therefore, G is a monotonically increasing 
function from G(0) = to G(x) > for all x G (0, a;*]. Therefore, the inequality fl30|) is 
proved and hence, the bound ffTTl) is verified. Note that ff29l) and fl30|) yield explicit bounds 
f[T2|) and f[T5]) . The proof of Theorem [T] is complete. 
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